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Abstract. The evolution with redshift of the temperature-luminosity relation of X-ray galaxy clusters is a key ingredient to 
break degeneracies in the interpretation of X-ray clusters redshift number counts. We therefore take advantage of the recent 
measurements of the temperature-luminosity relation of distant clusters observed with XMM-Newton and Chandra satellites to 
examine theoretical number counts expected for different available X-rays cluster samples, namely the RDCS, EMSS, SHARC, 
I60deg- and the MACS at redshift greater than 0.3. We derive these counts without any adjustment, using models previously 
normalized to the local (z ~ 0.05) temperature distribution function (TDF) and to the high-z (z ~ 0.33) TDF. We find that 
these models having Q.m in the range [0.85-1.] predict counts in remarkable agreement with the observed counts in the dilferent 
samples. We illustrate that this conclusion is weakly sensitive to the various ingredients of the modeling. Therefore number 
counts provide a robust evidence of an evolving population. A realistic flat low density model (Q.m = 0.3), normalized to the 
local abundance of clusters is found to overproduce cluster abundance at high redshift (above z ~ 0.5) by nearly an order 
of magnitude. This result is in conflict with the popular concordance model. The conflict could indicate a deviation from the 
expected scaling of the M - T relation with redshift. 
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1. Introduction 

The evolution of the temperature distribution function (TDF) 
of X-ray clusters is known to be a powerful cosmological test 
of the density parameter of the Universe (Oukbir & Blanchard 
1992). Indeed, the evolution with redshift of the mass function, 
on cluster scales, once normalized to present day, appears to 
be a robust cosmological test, with an exponential sensitivity 
to Q.M through the gravitational growth rate of perturbations 
(Blanchard & Bartlett 1998). The high sensitivity of this test 
has allowed its first direct application from a sample of ten (re- 
vised to 9) clusters at redshift 0.3 with measured temperature 
(Henry 1997). However, several authors have inferred some- 
what different values from this single sample (Eke et al. 1998; 
Viana & Liddle 1999; Blanchard et al. 2000, hereafter BOO). 
This might not be so surprising given the small size of the sam- 
ple as well as the limited range of redshift. An alternative ap- 
proach to track the evolution of the abundance of clusters is 
from their redshift distribution in X-ray samples for which the 



selection function is known. This procedure, reducing the time 
telescope investment, has been applied to the EMSS sample 
(Oukbir & Blanchard 1997; Reichart et al. 1999) and to the 
RDCS sample (Borgani et al. 1999, 2001) but the conclusion 
on Q.M depends on the possible evolution of the L - T rela- 
tion (Sadat et al. 1998, hereafter SB098; Borgani et al. 1999, 
Novicki et al. 2002). The various ingredients used in this mod- 
eling were not necessary accurately known, most noticeably 
the temperature mass {T - M) relation, and the luminosity- 
temperature (L - r) relation and its possible evolution with 
redshift. Therefore doubts have been raised up on the appli- 
cability of this test given these uncertainties (Colafrancesco et 
al. 1997; Rosati et al. 2002). 

The XMM-Q project was conducted in order to provide an 
accurate estimation of the possible evolution of the luminos- 
ity temperature relation at high redshift for clusters of medium 
luminosity which constitute the bulk of X-ray selected sam- 
ples, allowing to remove a major source of degeneracy in the 
determination of (Bartlett et al. 2001). To maintain a bet- 
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Table 1. Models and parameters used in the number counts cal- 
culations. 



ter control on the resulting L - T it is clearly preferable to use 
clusters homogeneously obtained from X-ray selected samples. 
Therefore, the choice of the SHARC surveys in the XMM- 
Q project provides several advantages: large angular coverage 
with the Bright SHARC (Nichol et al. 1999; Romer et al. 2000), 
deepness of the sample with the South SHARC (Burke et al. 
1997), stiU keeping the number of clusters to a reaUstic size for 
comprehensive X-ray investigations. 

The purpose of this letter is to examine the expected num- 
ber counts in comparison with the observed counts in two dif- 
ferent cosmological models, namely a concordance model and 
a high matter density flat universe, and to examine the ampli- 
tude of the major sources of uncertainties: the statistical dis- 
persion on the value of erg from the finite local sample, the 
systematic uncertainties in the M - T and the L - T relations 
at high redshift. In this modeling we take advantage of the ac- 
curate knowledge of the L - T relation provided by the first 
results from the XMM-f2 project (Lumb et al. 2003, hereafter 
L03) and from recent Chandra measurements of distant clusters 
(Vikhlinin et al. 2002). 

2. Modeling the temperature distribution function 
and number counts 

Although a fuU Ukehhood is possible in order to determine the 
best parameters and their final uncertainty taking into account 
all the possible sources of uncertainties, we find enlightening 
to illustrate the differences in the counts predicted in the two 
following specific models: the first model (hereafter model A) 
is the best flat model obtained as in BOO by fitting the local 
TDF and the high redshift TDF from the Henry (1997) sample, 
assuming a T CDM-like spectrum. The second model (here- 
after model B) we use is the so-called concordance model i.e. 
a flat A CDM model with Q.m - 0.3, in agreement with the 
recent WMAP results (Spergel et al. 2003). In this model, only 
the local TDF was then fitted. The parameters used are sum- 
marized in Table 1 . The samples we used for comparison are 
EMSS (Gioia et al. 1990; Henry et al. 1992), RDCS (Rosati et 
al. 1995, 1998), leOdeg^ (Vikhlinin et al. 1998, 2002), Bright 
SHARC (Nichol et al. 1999; Romer et al. 2000) and MACS 
(Ebeling et al. 2001) for which the selection criteria are be- 
lieved to be well known, generally given as sky coverage versus 
flux limit. These catalogs provide us with a sample of 274 (non 
necessarily independent) clusters with redshift between 0.3 and 
up to more than 1 and luminosities between 10^^ and 10*^ erg/s. 

In our modeling we follow a procedure close to the one 
used by Oukbir & Blanchard (1997), Reichart et al. (1999) and 
Borgani et al. (1999). In the first step, models are normalized 



using the local TDF, which request two fundamental ingredi- 
ents: the mass function and the M-T relation, assumed to fol- 
low a standard scaling law (Kaiser 1986). Intensive numerical 
simulations have allowed to provide accurate analytical fits to 
the mass function (Sheth, Mo & Tormen 2001, hereafter SMT; 
Jenkins et al. 2001; White 2002). Here we use the expression 
from SMT: 

— = J — C-— 1 + ^ exp (1) 

dm V n mdm \ (av^)i' ) \ 2 ) 

with a = 0.707, c = 0.3222 and p = 0.3 and v = 

2.1. The M-T relation 

It has been shown that the normalization of the M - T re- 
lation based on numerical simulations is significantly diff^er- 
ent from the normalization inferred from hydrostatic equation 
(Roussel et al. 2000). In this work we use two difl'erent normal- 
izations so as to cover the whole range of possibilities: we use 
on one side the calibration based on the numerical simulations 
of Bryan and Norman (1998, BN98 hereafter) and on the other 
side we use a M-T relation derived from the hydrostatic equa- 
tion (Markevitch 1998, M98 hereafter) which produces a lower 
normalization of the matter power spectrum o-g (Seljak 2002; 
Reiprich & Bohringer 2002). The M-T relation is written to 
be: 

T = Ti^i^MAiz, QM)/178)'/Xf (1 + ^) (2) 

where A(z, Qm) is the contrast density with respect to the den- 
sity of the Universe for virialized objects (hereafter h = 0.5). 
The subscript 15 means that masses are taken in unit of lO'^ so- 
lar masses. BN98 found T^ -3.8 keV, while M98 concluded 
to a higher normalization : ~ 6.5 keV. 

The mass function can then be normalized from the ob- 
served TDF (a possible dispersion of 15% is included as done 
in BOO). The availability of bright ROSAT clusters samples 
has allowed reasonably accurate estimations of the local TDF 
(M98; BOO; Pierpaoli et al. 2001; Ikebe et al. 2002), even if 
the agreement is partly due to the fact that the samples used in 
these previous works contain nearly the same clusters. 

In the present work, we use the local TDF based on an up- 
dated version of the local X-ray clusters sample used in BOO 
in which clusters were selected with fluxes fx > 2.2 X 10~'' 
erg/s/cm-^ in the ROSAT [0.1-2.4] keV and galactic coordinates 
l^iil > 20". This revised TDF is essentially identical to BOO and 
other recent measurements. 

2.2. Evolution of theL-T relation 

In a second step we use the L - T relation and its evolution 
to compute the luminosity function at diff'erent redshifts. Using 
the observed local L-T relation, models fitting the observed lo- 
cal TDF automatically provide a luminosity function matching 
the data (Oukbir, Bartlett & Blanchard, 1997). The local L-T 
relation Lboi - A(r/lkeV)* is estimated from the sample de- 
scribed above (Vauclair et al. 2003 in preparation). We obtained 
A = 0.0625 X 10^ erg/s assuming B = 3. We then estimated the 
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Fig. 1. XMM measurements of the evolution of clusters L-T 
relation expressed by the C(z) coefficient (Eq.|3}- Triangles are 
our XMM data and error bars are derived from the Icr error on 
the temperature measurement. Grey (yellow) area represents 
the Icr error on the C{z) fit and stars are the Chandra data. 

evolution of this relation using our XMM data of high-z clus- 
ters obtained in the XMM-O project (L03). Following SB098 
we compute for each cluster: 



C{z) 



AT'' Z),(Qm,z)2 
and fit these C(z) by a power-law evolution law: 

C(z)^(l+zf 



(3) 



(4) 



the best-fit parameter j8 - 0.65 + 0.21 was obtained by a stan- 
dard fitting (See figure ^ in good agreement with what 
was found previously (SB098) and fully consistent with the 
Chandra results (Vikhlinin et al. 2002). 



faint ROSAT samples number counts at redshift smaller than 
z = 0.3 for which clusters have too low luminosity compared 
to those used to measure the L-T evolution in this project. 
In addition at z < 0.3 clusters become significantly more ex- 
tended than the detection cell in the EMSS sample, and are not 
considered as reliable in MACS, a further reason for restricting 
our comparison to higher redshifts. 

2.4. Comparison with observations 

The final number counts for the different samples are presented 
in figure 121 for models A and B. As one can see the predicted 
counts in model A are in remarkable agreement with the ob- 
served number counts given the fact that no adjustment was 
performed, while the predicted number counts in the concor- 
dance model reveals a strong disagreement with the data. The 
overproduction of clusters at high redshift is impressively large, 
reaching a factor of nearly ten at redshift greater than 0.5. Two 
factors are responsible for this difference: the primary effect 
comes from the evolution in the number of clusters which is 
quite different in the two models and the second effect comes 
from the increase of the volume element in the A dominated 
universe. The uncertainties in the value of erg, in the L - T re- 
lation and in the M - T relation represent similar uncertainties 
in the predicted counts, in the range of 30% to 70 %, but are far 
below the difference between the two models which is between 
5 and 10. We also checked that changing the slope and the nor- 
malization of the local L-T does not affect our results. One 
may worry whether the selection procedure of the samples has 
introduced more bias than usually assumed. A systematic bias 
in the flux of the order of 2 to 3 for the different surveys would 
eliminate most of the discrepancy between the predictions of 
concordance model and the observed counts. Such a possibil- 
ity seems very unlikely: different flux inter-comparisons exist, 
including those of L03 between ROSAT SHARC and XMM 
fluxes and those of Vikhlinin et al. (2003) between ROSAT and 
Chandra fluxes, and do not provide any evidence for such a 
large bias. 



2.3. Number counts 

In order to compute number counts, one can notice that the 
observations actually provide z and fx (rather than the actual L^ 
and T-t). For a flux limited sample with a sky coverage area Q. 
and a flux limit /( one has therefore to compute the following: 

A^O /,, z, Az) = Q j^^ f (L. > 47TDffx)dz 

= £lJ^^N(>T(z))dV(z) 

= " i-A? ImZ ^')dMdV{z) (5) 

where T(z) is the temperature threshold corresponding to the 
flux fx as given by the observations, being therefore indepen- 
dent of the cosmological model. We then compute the theoret- 
ical number counts in redshift bins Az = 0. 1 for the various 
clusters samples mentioned above accordingly to their respec- 
tive selection function. In order to compare the same cluster 
population at different redshifts, we have excluded from the 



3. Discussion 

As we have seen a model which is normalized to the local 
as well as to the high redshift TDF, reproduces impressively 
well the redshift distribution of all the surveys we have investi- 
gated, without any adjustment and with little uncertainties aris- 
ing from the modeling. In contrast, we have checked that relax- 
ing the request of matching the local abundance a concordance 
model fitting the redshift distribution leads to a local abundance 
of X-ray clusters A^(> 4keV,z - 0.05) nearly 5 times lower 
than actually observed. This is a strong indication that existing 
samples of clusters (namely the Henry sample, the RDCS, the 
EMSS, the Bright SHARC, the 160deg2 and the MACS sam- 
ple) draw the same picture, consistently pointing out toward the 
fact that the cluster abundance is significantly evolving with 
redshift, perhaps closing a long term controversy on this ques- 
tion (Henry et al. 1992). Furthermore, in standard hierarchical 
picture of structure formation such evolution points toward a 
high matter density universe with Qm in the range [0.85-1.], 
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the precise value depending on the M - T normahzation. This 
conclusion is clearly conflicting with the currently popular con- 
cordance model. However, it should be emphasized that this 
is entirely consistent with all previous existing analyzes per- 
formed on the redshift distribution of X-ray selected samples 
of clusters performed with the same methodology: our conclu- 
sion corroborates previous analyzes of the EMSS clusters red- 
shift distribution: SB098 found Qm - 0.85 + .2, Reichart et 
al. (1999) found Q.m = 0.96 ± .3 as weU as the RDCS red- 
shift distribution for models normalized to the present-day X- 
ray clusters abundance. Indeed, from Eq. (14) in Borgani et al. 
(1999), Qm = (j6 + 2)/3 ± 1/3 we derive Qm ~ 0.88 ± .34 from 
our yS ~ 0.65 + .21. Note that high Qm models are also consis- 
tent with WMAP (Blanchard et al. 2003). A possible loophole 
could be a large systematic bias flux, but we have argued that 
this seems very unlikely. An other possibility would be that the 
scaling in the redshift of the M - T relation (Eq. 1) is com- 
pletely wrong, violating the basic scaling scheme. Voit (2000) 
has investigated such a possibility, but concluded to a moder- 
ate effect. In Fig. 2, we have plotted the predicted counts in a 
concordance model, assuming 

T = ri5(OMA(z, QM)/178)'/Xf (6) 

instead of Eq. 2. As one can see, such a modification reestab- 
lishes agreement of the concordance model with observations. 
It is well known that the L - T relation cannot be explained 
from simple scahng arguments. One may therefore argue that 
the redshift evolution of the M - T relation may suffer from 
more dramatic effect than usually assumed, although - to our 
knowledge - such a possibility has never been advocated and it 
is probably not obvious to find physical motivation leading to 
gas thermal energy in distant clusters (z ~ 1) to be reduced by 
a factor of two compared to clusters in the local universe. We 
therefore conclude that the redshift distributions of present-day 
available X-ray clusters surveys, as well as the recent results on 
the L-T relation of high redshift clusters, favor a high matter 
density universe unless the standard paradigm on clusters gas 
physics has to be deeply revised. 
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Fig. 2. Theoretical number counts in bins of redshift (Az = 0.1) for the different surveys: RDCS, EMSS, SHARC, MACS and 
160deg^. Observed numbers are triangles with 95% confidence interval on the density assuming poissonian statistics (arrows 
are 95% upper limits). For the 160deg^ survey we have also examined the counts for the brightest part (160deg^(high flux) 
corresponding to fluxes /r > 2 x 10"'-' erg/s/cm^), in order to show that the excess seen at z ~ 0.5 is not due to luminous clusters. 
The upper (red) curves are the predictions in the concordance model (model B). The lower (blue) curves are for model A (see 
table 1.). Diff'erent M - T are figured: the dashed lines correspond to Tig - 6.5keV (M98 M - T) while the continuous lines 
are for Tig = 4.keV, corresponding to virial mass nearly twice larger, close to BN98 normalization. The dark grey (green) area 
is the uncertainty range from our estimates on the uncertainty in the evolution of the L - T relation while the light grey area 
is the uncertainty on number counts due to the dispersion on erg. The 3-dotted-dashed lines show the predicted counts in the 
concordance model using M - T relation from Eq. 6 violating the standard scaling with redshift. 



